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Abstract 



Q ■ The spin-1/2 Ising-Heisenberg model on diamond-like decorated Bethe lattices is exactly solved 

with the help of decoration-iteration transformation and exact recursion relations. It is shown that 

the model under investigation exhibits reentrant phase transitions whenever a sufficiently high 

. ' coordination number of the underlying Bethe lattice is considered. 
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I. INTRODUCTION 

The ferromagnetic spin-1/2 Ising-Heisenberg model on planar lattices currently attracts 
a great deal of attention as it may exhibit a variety of interesting phenomena such as a 
quantum critical point [l|, l2|, peculiar spontaneous antiferromagnetic long-range order [l(] 
or several unusual spin-liquid ground states 2|. In this work, the spin-1/2 Ising-Heisenberg 
model on diamond-like decorated Bethe lattices will be treated exactly with the aim to study 
its ground-state and finite-temperature phase diagrams. 

II. ISING-HEISENBERG MODEL AND ITS EXACT SOLUTION 

Consider the spin-1/2 Ising-Heisenberg model on the diamond-like decorated Bethe lat- 
tices schematically illustrated on the left-hand-side of Fig. [H In this figure, the full circles 
denote lattice positions of the Ising spins fi, the empty circles label lattice positions of the 
Heisenberg spins S and the parameter q is the coordination number of the underlying Bethe 
lattice. For further convenience, the total Hamiltonian can be defined as a sum of the bond 
Hamiltonians, i.e. 'H = 'Ylik'^ki where each bond Hamiltonian l-i^. involves all the interaction 
terms belonging to the kth. diamond unit (see Fig. [1]) 

T-ik = — J[^{Sj^iSj^2 + Skl^k2) + ^klSk2\ 

-Ji{si, + si,Wk, + iii,) (1) 

In above, S^ (a = x,y,z) and fij label spatial components of the spin-1/2 operator, the 
former term J accounts for the XXZ interaction between the nearest-neighbour Heisenberg 
spins, A is a spatial anisotropy in this interaction and the latter term Jj accounts for the 
Ising-type interaction between the nearest-neighbour Ising and Heisenberg spins, respec- 
tively. 

A crucial step of our approach lies an evaluation of the partition function. Due to 
a commutability of different bond Hamiltonians, the partition function can be partially 
factorized to a product of bond partition functions 

Nq/2 Nq/2 

Zmu = 5Z n Tr-fc expi-pn,) = 5^ n ^fc, (2) 

{fJ-i} k=l {^J fc=l 

where /3 = l/{k-BT), k^ is Boltzmann's constant, T denotes absolute temperature, A^ is the 
total number of the Ising spins, the symbol Tryt stands for a trace over degrees of freedom of 
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FIG. 1: A schematic representation of the decoration- iteration transformation (DIT) between the 
spin-1/2 Ising-Heisenberg model on the diamond-hke decorated Bethe lattice (figure on the left) 
and the spin-1/2 Ising model on the simple Bethe lattice with the coordination number q = 3 
(figure on the right). The full (empty) circles denote lattice positions of the Ising (Heisenberg) 
spins and the ellipse demarcates all interaction terms belonging to the fcth bond Hamiltonian ([!]). 

the A;th Heisenberg spin pair and the summation X]|m I ^^ns over all available configurations 
of the Ising spins. A calculation of the bond partition function Z^ can easily be accomplished 
through a direct diagonalization of the bond Hamiltonian ([1]). Subsequently, the gained 
expression for Z^ will depend just on two outer Ising spins and can be replaced with the 
generalized decoration-iteration transformation [2| 

Mf^li, /^fe2) = 2 exp(/3 J/4) cosh[/3 Ji(/i^i + fxl^)] + 

2exp(-/3J/4) cosh(/3JA/2) = Aexp{/3Rfxl^fxl2)- (3) 

The mapping parameters A and R can directly be obtained from the decoration-iteration 
transformation ([3]), which holds quite generally if and only if 

A=[Z,{l/2,l/2)Z,{l/2,-l/2)f\ (4) 

/3R = 2 \n[Zk{l/2, 1/2)] - 2 ln[Z,(l/2, -1/2)]. (5) 

Substituting the algebraic transformation ([3]) into the formula (|2]) one gets a precise map- 
ping relationship between partition functions of the spin-1/2 Ising-Heisenberg model on the 
diamond-like decorated Bethe lattice and the corresponding spin-1/2 Ising model on the 
simple Bethe lattice (see Fig. [1]) 

Zihm(/3, J, A, Ji) = A^''/2Zim(/3, R). (6) 

Note that the partition function of the spin-1/2 Ising model on the Bethe lattice can exactly 
be calculated by making use of exact recursion relations [sl and hence, our exact calculation 



is formally completed. Using the mapping equivalence ([6]) between both partition functions, 
the sublattice magnetization of the Ising spins directly equals to the single-site magnetization 
of the corresponding spin-1/2 Ising model on the simple Bethe lattice 

^i = (Afc)iHM = if^DiM = miMi(3R). (7) 

Above, the symbols (■ ■ ■)ihm and (■ ■ ■)im mark canonical ensemble average performed within 
two models connected through the mapping relation (Q. Note furthermore that the corre- 
sponding exact result for the magnetization ttijm can be also found with the help of recursion 

n 

relations ^] . The sublattice magnetization of the Heisenberg spins follows from the formula 

^z_/^z^ miM sinh(/3 Ji) 

cosh{(3 Ji) + exp(-^) coshf ^j 

III. RESULTS AND DISCUSSION 

First, let us take a closer look at the ground state of the spin-1/2 Ising-Heisenberg model 
with both ferromagnetic interactions (J > 0, Ji > 0). The classical ferromagnetic phase with 
all Ising as well as Heisenberg spins aligned parallel constitutes the ground state whenever 
the exchange anisotropy A < Ac = 1 + 2Ji/J. In the opposite case A > Ac, the remarkable 
spin-liquid phase constitutes the ground state and this disordered phase can be character- 
ized through a complete randomness of the Ising spins and the quantum entanglement of 
Heisenberg spin pairs. An origin of this unusual disordered phase lies in the geometric frus- 
tration, which is triggered by a competition between the easy-axis Ising interaction and the 
easy-plane XXZ Heisenberg interaction 2J]. 

Fig. [2] shows the critical temperature as a function of the exchange anisotropy A for one 
particular ratio J/Ji = 1 and several values of the coordination number q. As one can 
see, the critical temperature vs. exchange anisotropy dependence basically changes with the 
coordination number q of the underlying Bethe lattice. The critical lines are approaching the 
zero-temperature phase boundary between the classical ferromagnetic phase and disordered 
phase at Ac with a negative (positive) slope for the diamond-like decorated Bethe lattices 
with the coordination number g < 4 (g > 4) and with an infinite gradient for the particular 
case with the coordination number g = 4. Owing to this fact, reentrant phase transitions 
can be observed for A > Ac only when assuming the diamond-like decorated Bethe lattices 
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FIG. 2: Critical temperature as a function of the exchange anisotropy A for one selected ratio 
J/Ji = 1 and several values of the coordination number q. 
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FIG. 3: Thermal variations of both spontaneous sublattice magnetizations for J/Ji = 1, g = 6 and 
two different values of the exchange anisotropy A. 

with a sufficiently high coordination number g > 4. Thermal variations of both sublattice 
magnetizations, which are depicted in Fig. [3l provide an independent confirmation of the 
aforedescribed reentrance. 

In conclusion, it is worthy to note that the rigorous procedure developed on the grounds of 
decoration-iteration transformation and exact recursion relations can be readily adapted to 
treat the investigated model in a presence of external magnetic field, which will be examined 
in detail in our forthcoming work. 
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